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ABSTRACT 

The  Final  Report  details  the  research  results  for  the  full 
period  from  August  2001  to  July  2002.  It  includes  a  brief 
overview  on  chaotic  electronic  oscillators  operating  in  the 
RF  and  the  microwave  ranges.  Special  attention  is  paid  to 
the  delay  line  based  chaotic  oscillators.  The  original 
research  contributes  to  the  numerical  and  experimental 
investigation  of  chaotic  and  hyper-chaotic  oscillators  for  the 
very  high  frequency  (VHF:  30  to  300  MHz)  and  the 
ultrahigh  frequency  (UHF:  300  to  1000  MHz)  ranges.  The 
design  of  three  different  versions  of  a  delay  line  based 
oscillator,  the  development,  linear  analysis  and  numerical 
simulation  of  the  appropriate  mathematical  models, 
eventually  building  and  experimental  investigation  of  the 
hardware  prototypes  are  described. 


1.  INTRODUCTION 

The  first  electronic  chaos  oscillators  began  to  emerge  in  the  late  70s  and  early  80s. 
Electronic  circuits  exhibiting  chaotic  behavior  were  rather  cheap  and  relatively  easy  to  build. 
They  were  employed  as  convenient  tools  for  modeling  and  studying  chaotic  phenomena,  e.g. 
universal  routes  to  chaos,  control  of  chaos,  etc.  Meanwhile,  electronic  chaos  oscillators  can  be 
readily  exploited  as  high  power  and  high  efficiency  dynamic  noise  generators.  Another  burst 
of  the  interest  in  chaos  generators  occurred  in  the  early  90s  when  the  American  scientists  of 
Naval  Research  Lab.  in  Washington,  D.C.,  also  of  the  Massachusetts  Institute  of  Technology 
suggested  an  intriguing  idea  to  make  use  of  chaos  in  secure  communications.  However,  most 
of  the  electronic  chaos  oscillators  developed  from  the  80s  were  designed  and  investigated  at 
relatively  low  frequencies,  in  the  kilohertz  range.  Only  few  examples  of  dynamic  chaos 
generators  were  developed  to  operate  in  the  microwave  range,  e.g.  [1-8].  Meanwhile,  modem 
information  and  military  technologies,  including  computer,  radar  and  communication  systems 
need  very  high  frequency  (VHF:  30  to  300  MHz)  and  ultrahigh  frequency  (UHF:  300  to  1000 
MHz)  chaos  generators  to  be  developed  and  used  as  high  power  controllable  noise  sources. 
For  example,  in  1996  the  Mathematical  and  Computer  Sciences  Division,  the  Electronics 
Division  and  the  Physics  Division  of  the  U.S.  Army  Research  Office  (ARO)  held  an 
Advanced  Concept  Workshop  “ Communicating  by  Chaos:  Digital  Signal  Generation  by 
Simple  Nonlinear  Devices ”.  One  of  the  conclusions  of  the  Strategic  Assessment  Report  says 
that  “...chaotic  oscillators  for  use  as  a  transmitter  at  high  megahertz  to  gigahertz  frequencies 
are  urgently  needed  and  a  better  understanding  of  the  physics  of  nonlinear  and  active 
components  is  required”  [9]. 
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Dynamical  systems  with  time  delays  exhibit  very  rich  dynamics,  including  high¬ 
dimensional  hyper-chaotic  behavior  characterized  by  multiple  positive  Lyapunov  exponents, 
broadband  continuous  spectra,  thus  can  be  exploited  as  noise  like  signal  generators.  These 
systems  are  commonly  presented  by  the  delay  differential  equations  (DDE): 


^!l  =  F[x(t),x(t-  T)]. 


(1.1) 


Here  x(t)  is  a  dynamical  variable,  t  is  time,  ris  time  delay,  and  F[...]  is  a  nonlinear  function. 
A  specific  example  of  this  type  of  DDE  provides  the  Mackey-Glass  (MG)  mathematical 
model  [10]: 


dx(t)  _  2x(t  —  T ) 

~dT~  l  +  xl0(t-T)~X(t 


(1.2) 


The  first  right-hand  term  in  the  Eqn.(2)  represents  a  nonlinear  function  with  delay.  The  MG 
model  has  been  studied  numerically  in  details  by  Fanner  [11]. 

To  our  best  knowledge  an  electronic  circuit  imitating  the  MG  system  has  been  first 
described  by  Dmitriev  et  al.  [12].  An  acoustic-electric  device  has  been  employed  for  the  delay 
element.  Later  an  analog  computer  design  has  been  used  by  Losson  et  al.  [13]  to  simulate  the 
dynamics  of  a  related  model  -  the  piecewise  constant  version  of  the  MG  model  containing  a 
rectangular  nonlinear  function.  The  delay  unit  in  their  analog  computer  has  been  implemented 
by  means  of  a  sample-and-hold  (S&H)  chip,  or  the  so-called  “bucket  brigade  device”  used 
along  with  a  clock  generator.  Actually  the  computer  is  a  hybrid  one  containing  both  time 
continuous  and  time  discrete  units.  Later  a  DDE  with  a  tent  type  nonlinear  function  as  well  as 
a  rectangular  one  have  been  considered  by  Schwarz  and  Mogel  [14].  In  addition,  the  authors 
of  this  paper  suggested  to  replace  the  S&H  device  by  a  transmission  line  (e.g.  a  coaxial 
cable),  thus  yielding  a  fully  time  continuous  system.  Another  version  of  a  time  continuous 
delay  system  has  been  proposed  by  Namajunas  et  al.  [15].  Two  complementary  junction  field 
effect  transistors  (JFET)  coupled  in  a  special  way  [16]  were  used  to  compose  the  nonlinear 
device  while  the  delay  unit  was  implemented  by  means  of  a  high-order  L-C-L  low-pass  filter. 
The  nonlinear  device  composed  of  the  complementary  pair  of  JFET  is  described  by  a  function 
very  similar  in  fonn  to  that  introduced  in  the  MG  model.  Therefore,  this  analog  circuit  [15] 
was  employed  to  investigate  various  features  of  the  hyper-chaotic  MG  system,  including 
control  of  chaos/hyperchaos  [17,18]  and  synchronization  of  chaos/hyperchaos  [19-21].  In 
addition,  Pyragas  has  considered  the  possibility  to  apply  the  MG  hyper-chaotic  oscillators  to 
secure  communications  [22]. 

Most  experiments  with  the  MG  and  related  type  electronic  circuits  have  been  performed 
at  relatively  low  (kilohertz)  frequencies  [12,13,15,17-21].  High  frequency  delay  line  based 
chaotic  oscillator  has  been  considered  by  Mogel  et  al.  [23].  However,  the  paper  includes  only 
a  short  note  that  a  practical  prototype  generates  noise  like  signals  in  a  spectral  range  from  DC 
to  over  100  MHz.  Neither  details  about  the  nonlinear  units  nor  experimental  power  spectra  are 
presented  in  [23]. 
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2.  DELAY  LINE  OSCILLATOR 

The  diagram  of  a  delay  line  oscillator  is  sketched  in  Fig.  1.  The  oscillator  has  a  ring 
structure  and  contains  five  main  units:  (1)  a  nonlinear  element  NE,  (2)  an  RC  low-pass  filter, 
(3)  an  amplifier  AMP,  (4)  a  Delay  line,  and  (5)  a  buffer  BUF.  The  latter  is  an  auxiliary  unit 
and  may  be  skipped  in  some  specific  designs.  The  Rl  is  an  external  load  (commonly  50  Q). 


kx(  6-  r) 


Fig.  1.  Block  diagram  of  a  delay  line  based  oscillator. 


2.1.  Nonlinear  element 


A  simplified  circuit  diagram  of  a  single  transistor  Q1  based  nonlinear  element  and  its 
input-output  function  Vc  vs.  VB  are  shown  in  Fig.2. 


vi 


Fig.  2.  Circuit  diagram  of  the  nonlinear  element  and  its  input-output  function  for  b>a. 


The  input  (base)  voltage  Vb  is  a  full  signal  (dc  bias  Vbo  +  ac  signal  V).  The  dc  component  of 
the  output  (collector)  voltage  appearing  due  to  the  Vi  is  removed  for  simplicity.  Thus  Vc 
corresponds  to  an  ac  component  only.  Until  Vg<U*  (V*  is  the  breakpoint  voltage  in  the 
current-voltage  characteristic  of  the  base-emitter  junction;  V*~0.7V  for  silicon  transistors)  the 
transistor  Q1  is  in  the  OFF-state  and  the  output  signal  is  due  to  the  by-pass  resistor  Rl  only: 


yc=. 


R, 


R\  +Ri 


VB  =  aV 


B  • 


When  VB>V*  the  Q1  is  in  the  ON-state  and 
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yc  = 


— — — VB  V  -V*)  =  aVB  -b(VB  - V *). 

Rx+R2  R3  +  r 


Here  R1WR2  is  the  resistance  of  Rl  and  R2  connected  in  parallel,  and  r  is  the  small  signal 
resistance  of  the  base-emitter  junction  of  the  Q1  in  the  ON  state  near  the  operating  point  O.  If 
b>a  the  falling  segment  of  the  input-output  (transfer)  function  does  appear.  It  is  convenient  to 
normalize  all  the  voltages  of  interest  to  the  breakpoint  voltage  V *:  Vc/V*=F,  Vb/V*=v, 
Vbo/V*=c,  and  V/V*=u.  Then  the  transfer  function  is  given  by: 


F(v)  = 


av, 

(1 a  -  b)v  +  b, 


v<l, 

V  >  1 . 


(2.1) 


Since  v=c+u  the  transfer  function  can  be  presented  in  the  following  form: 


FA(n) 


j  au  +  ac,  c  +  u  <  1, 

[(a  —  b)u  +  ac  —  b(c  —  1),  c  +  u>  1. 


(2.1a) 


Now  the  operating  point  O  appears  on  the  vertical  axis  of  the  plane  [F,  u\.  In  addition,  the 
operating  point  O  (at  c>l)  can  be  moved  down  to  the  origin  by  eliminating  the  constant  term 
ac-b(c-  \ ): 


F(u)  = 


\au  +  b(c  - 1), 
|  ( a  -  b)u, 


c  +  u<  1, 
c  +  u  >  1 . 


(2.1b) 


2.2.  RC  low-pass  filter 

The  lst-order  RC  low-pass  filter  shown  in  Fig.3. 


R 

Ui„  * - 1  | — • - m  Uci 


Cl 


Fig.  3.  Circuit  diagram  of  the  1st-order  RC  filter. 


Dynamics  of  the  filter  is  described  by  the  linear  ordinary  differential  equation  (ODE): 

c  dUci  _  Uin  -Uci 


dt 


R 


(2.2) 


Here  Uci  is  the  voltage  across  the  capacitor  Cl.  Note,  that  in  a  specific  circuit  the  output 
resistance  of  the  previous  unit,  for  example  R\\\R2  of  the  nonlinear  element  may  play  the  role 
of  R.  By  introducing  the  dimensionless  variables 


x  = 


Uci_ 
V *  ’ 


x 


(lx 

~d6 


and  considering  the  input  voltage  Um  to  be  just  the  output  voltage  the  nonlinear  element,  i.e. 
U in~V c  Eqn.(2.2)  can  be  given  in  the  following  fonn: 

x  =  F(u)-x.  (2.3) 


2.3.  Amplifier 

This  unit  is  a  linear  wideband  amplifier  with  the  gain  k: 


U out  kU in  . 


Since  the  input  of  the  amplifier  Uu,  is  just  the  output  of  the  RC  filter,  Uin=Uci  the  normalized 
output  X=U0JV*  is 


X  =  kx  .  (2.4) 

The  amplifier  should  be  a  non-inverting  one,  therefore  it  should  contain  an  even  number  of 
stages  (minimum  two). 


2.4.  Delay  line 

Any  transmission  line,  say  a  segment  of  50  Ohm  coaxial  transmission  cable  can  be  used 
for  the  delay  line  (time  delay  Tdei  is  about  5  ns  per  meter).  The  end  of  the  transmission  line 
should  be  loaded  with  a  matched  resistor  to  avoid  reflected  signals.  Taking  into  account 
formula  (2.4)  the  normalized  delayed  signal  is  given  by 


Xdel=X(0-T)  =  kx(d-Z). 


(2.5) 


Here  the  parameter  ris  a  reduced  time  delay,  T  =Tdei/RC \ . 


2.5.  Buffer 

The  buffer  unit  is  a  simply  one-  or  two-stage  emitter  follower.  Due  to  its  high  input 
impedance  it  “isolates”  the  Delay  line  from  the  nonlinear  unit  NE,  thus  serves  to  avoid 
undesirable  loading  of  the  Delay  line.  In  addition,  because  of  its  low  output  impedance  the 
external  load  Rl  does  not  effect  the  dynamics  of  the  oscillator. 
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3.  Delay  differential  equations 

3.1.  Case  1:  Single  oscillator 

By  substituting  the  indefinite  variable  u  in  Eqn.(2.1b)  and  in  the  ODE  (2.3)  with  Xdei 
from  relation  (2.5)  we  obtained  a  DDE  describing  the  closed-loop  (ring)  circuit  in  Fig.l: 


x  =  F[x(d  -t)]  -  x. 

Here 

jakx(0  —  t)  +  b(c  - 1),  c  +  kx(6  -t)  <  1, 
[  (a  -b)kx{6  -r),  c  +  kx{6  —  z)  >  1. 


The  nonlinear  function  F  is  depicted  for  illustration  in  Fig.4. 

F(x) 


(3.1) 

(3.2) 


Fig.  4.  Nonlinear  transfer  function  F{x).  a= 0.5,  b=  1.8,  c=l.l,  k= 3. 

At  T=0  Eqn.(3.1)  for  ak>  1  and  b>a  has  two  different  steady  state  solutions: 

__b(c-\) 

v0i  -  -  —  ana  x02  -  o . 

ak  - 1 

The  xo\  is  an  unstable  one,  while  the  vo2  at  the  origin  is  a  stable  one.  For  non-zero  delays 
Eqn.(3.1)  has  been  solved  numerically  using  the  Turbo  Pascal  software  and  the  Runge-Kutta 
integration  method.  At  small  delays  t<0.49  (for  the  given  parameters)  the  latter  remains 
stable.  At  z>0.49  it  turns  to  be  unstable  yielding  periodic  oscillations.  In  the  interval 
1.23<t<1.45  the  system  undergoes  period  doubling  bifurcations.  With  further  increase  of  Tit 
exhibits  chaotic  behavior  as  illustrated  in  Fig.  5. 


Fig.  5.  Typical  chaotic  waveform  from  Eqn.(3.1).  a=0.5,  b=  1.8,  c=l.l,  k= 3,  t^8. 
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3.2.  Case  2:  Two  “slow”  oscillators  in  a  ring 

Let  us  consider  more  complicated  case  when  two  oscillators,  modules  X  and  Y  each 
described  by  Eqn.(3.1)  are  coupled  in  a  ring  as  shown  in  Fig. 6.  “Slow”  oscillators  mean  that 
both  of  them  possess  the  RC  low-pass  filters  in  their  circuits,  that  is  the  C 1  in  the  X-module 
as  well  as  a  corresponding  capacitor  C2  in  the  Y-module  have  finite  non-zero  values. 


X 

Y 

Fig.  6.  A  system  of  two  coupled  delay  line  oscillators. 

The  coupled  system  is  described  by  a  set  of  two  coupled  DDE: 

x  =  F[y(0-z  i)]-*, 

£y  =  F[x(d-T2)\-y, 


(3.3) 


with  the  same  type  of  the  nonlinear  functions  F(x)  and  F(y)  as  given  by  (3.2).  In  Eqn.(3.3) 
e=Cj]C\,  where  C2  is  the  capacitance  of  the  RC  filter  in  the  Y-module  (in  general,  C&C\). 
Numerical  analysis  shows  that  two  delays  X\  and  To  can  be  replaced  by  their  sum  T=Ti+To, 
consequently  the  set  of  two  DDE  is  substituted  by  a  simpler  set  of  one  DDE  and  one  ODE: 


x  =  F[y(d-z)\-x, 

or 

ey  =  F(x)-y, 


x  =  F(y)-x, 

£y  =  F[x(6-z)]-y, 


(3.3a) 


In  (3.3a)  both  equations,  the  DDE  and  the  ODE  are  nonlinear  ones.  Making  use  of  the  fact 
that  in  the  ring  structure  any  element  (at  least  mathematically)  can  be  moved  to  any  place  in 
the  loop  the  set  (3.3)  can  be  further  simplified  and  presented  in  the  fonn  of  a  nonlinear  DDE 
and  a  linear  ODE: 


x  =  ®[y(d- t)]-x, 

or 

ey  =  x-y. 


x=y-  x, 

ey  =  <J>[x(0-T)]-y, 


(3.3b) 


Here  <t>=F {F[u(9-t)]} ,  i.e.  is  a  more  complicated  (double)  transfer  function  of  u: 

\aF  +  b(c  —  1),  c  +  F<l, 
d>  =  < 

I  ( a-b)F ,  c  +  F  >  1. 


(3.4) 


The  graphical  view  of  <f>(.r)  obtained  numerically  is  illustrated  in  Fig. 7.  It  is  a  four  segment 
(“double-camel”)  piecewise  linear  function.  Also  an  analytic  form  of  <J>(.r)  can  be  derived 
from  (3.2)  and  (3.4)  by  means  of  simple  algebra: 

a2k2x  +  b(c  -  l)(ak  + 1),  x<xl, 

^  _  a  (a  —b)kx  +  (a  —b)kb(c  - 1),  x\<x<x2, 

( a-b)2k2x ,  x2^x<x-i, 

a(a  -  b)k  x  +  b(c  - 1),  x>x2, 


(3.4a) 
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<*>(*) 


Fig.  7.  Nonlinear  transfer  function  <F(.x;).  a=0.5,  b=  1.3,  c=l.l,  k= 3.2. 


where  the  three  break-points  are  the  following: 

_M  +  l(c_1)^  x  x= - I — -(c-1). 

nk2  k  (, a-b)k 2 


For  the  given  parameter  values  Vi,  a'2<0  and  a'3>0  (<:/</?)  in  agreement  with  the  graph  in  Fig. 7. 


3.3.  Case  3:  “Slow”  and  “fast”  oscillators  in  a  ring 

In  the  limit  when  £  =0  (C2=0),  that  is  the  second  oscillator  is  “fast”  we  come  back  to  a 
single  DDE  similar  to  Eqn.(3. 1)  for  the  Case  1  but  with  a  different  nonlinear  function  : 


x  =  d>[.\(0  —  t)]  -  x.  (3.5) 

At  1^=0  Eqn.(3.5),  with  <£>(v)  described  by  (3.4  or  3.4a)  has  four  different  steady  state 
solutions: 


*oi  - 


in,  !) 

ak  - 1 


(a  -b)kb(c  - 1) 
a(a  -b)k2  - 1 


v03  =  0 ,  and 


*04  -  - 


b(c- 1) 
a(a  -  b)k 2  - 1 


The  xoi  and  the  V03  are  exactly  the  same  as  the  xo\  and  the  .v'02  in  the  Case  1  of  the  two  segment 
function  F(x).  However,  in  contrast  to  the  Case  1  they  are  both  unstable  solutions,  including 
V03  at  the  origin,  because  they  appear  on  the  raising  segments  of  cF(a).  Meanwhile  another  pair 
of  the  steady  state  solutions,  V02  and  *04  that  appear  on  the  falling  segments  of  the  transfer 
function  are  stable  ones.  Similarly  to  the  Case  1  they  become  unstable  for  non-zero  delays. 

The  steady  state  analysis  perfonned  here  at  t=0  for  the  Case  3  can  be  useful  for  the 
Case  2  as  well,  since  four  values  of  the  steady  state  solutions,  xq,  (i=l,..,  4)  are  identical  ones 
in  both  cases.  While  another  four  values,  yoi  are  obtained  immediately:  yoi=F(xoi)  for  Eqn.(3.3 
or  3.3a)  and  yoi=^oi  for  Eqn.(3.3b). 
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3.4.  Numerical  results 

In  this  section  the  qualitative  and  quantitative  numerical  results  for  all  the  three 
discussed  cases  are  presented.  The  results  include  phase  portraits,  Poincare  sections, 
Lyapunov  exponents,  Kolmogorov  entropy  and  power  spectra. 


Fig.  8.  Phase  portraits  (left)  and  Poincare  sections  (right)  in  three  different  cases: 
(a)  case  1  from  Eqn.(3.1).  a=0.5,  b=  1.8,  c=l.l,  k= 3,  t=8. 

(b)  case  2  from  Eqn.(3.3a).  a= 0.5,  b=  1.8,  c=l.l,  k=  3,  x=  8,  e=l. 

(c)  case  3  from  Eqn.(3.5).  a= 0.5,  b=  1.3,  c=l.l,  k= 3.2,  t=8,  (e=0). 

The  Poincare  sections  are  taken  at  the  local  maxima  of  the  x(t),  that  is  at  dx/dt=  0. 
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One  can  see  that  the  phase  portraits  in  Case  2  and  Case  3  are  more  complicated  than  in 
Case  1,  that  is  similar,  as  expected,  to  the  phase  portrait  of  the  MG  system  [11,15],  An 
interesting  fact,  that  the  Poincare  sections  in  Case  1  and  Case  3  are  linear  (one-dimensional) 
ones.  Moreover,  they  do  coincide  with  the  transfer  functions  F(x )  in  Fig.4  and  with  <E(.i')  in 
Fig. 7,  respectively.  However,  this  is  not  an  unexpected  feature.  It  comes  in  a  straightforward 
way  from  Eqn.(3.1)  and  Eqn.(3.5).  Indeed,  if  we  set  the  Poincare  condition  dx/dt= 0,  then 
x=F[x(0- t)]  and  v=<I>[.r(0-T)]  in  the  Case  1  and  Case  3,  respectively.  Meanwhile  the  dots  of 
the  Poincare  section  in  the  Case  2  is  scattered  in  the  plane,  that  is  an  essentially  two- 
dimensional  one. 

In  addition,  all  the  three  cases  have  been  characterized  quantitatively  by  the  Lyapunov 
exponents  (LE)  (Fig. 9).  In  Case  1  the  system  becomes  chaotic  (has  one  positive  LE)  at  t=1.5 
as  discussed  above  on  page  9.  With  the  increase  of  z  more  and  more  positive  LE  do  appear, 
that  is  the  system  becomes  a  hyper-chaotic  one.  Similar  behavior  is  observed  in  Case  2  and 
Case  3.  The  difference  is  in  the  fact,  that  chaos  emerges  at  lower  delays  (t=1)  and  the  values 
of  the  LE  are  nearly  twice  larger. 


Fig.  9.  Non-negative  Lyapunov  exponents  A,  versus  the  delay  Tin  three  different  cases: 
(a)  case  1 ,  (b)  case  2,  (c)  case  3.  Parameters  are  the  same  as  in  Fig. 8. 
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The  results  of  Fig.9  are  summarized  in  Fig.  10,  where  the  number  of  positive  LEs  n+  and 
the  Kolmogor  entropy,  i.e.  the  sum  of  the  positive  LE:  H=L2?  are  presented. 


T 


Fig.  10.  Number  of  positive  Lyapunov  exponents  (left)  and  the  Kolmogorov  entropy  (right). 
The  letters  a,  b,  and  c  in  the  legends  correspond  to  the  cases  of  Fig.9. 

Finally,  we  present  in  Fig.  1 1  the  power  spectra  simulated  from  the  variable  x(t).  One 
can  see  that  the  most  smooth  spectrum  is  in  the  Case  2. 


Fig.  11.  Simulated  power  spectra  for  three  different  cases: 

(a)  case  1 ,  (b)  case  2,  (c)  case  3.  Parameters  are  the  same  as  in  Fig. 8. 
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4.  Hardware  prototypes  and  experimental  results 

We  have  developed  three  practical  modifications  of  the  delay  line  based  chaotic 
oscillator:  (?)  a  “DC”  one  with  the  bandwidth  of  its  units  from  dc  to  several  gigahertz,  (ii)  a 
broadband  “AC”  one  with  the  bandwidth  of  its  units  from  several  megahertz  to  several 
gigahertz,  and  (iii)  an  impoved  version  of  the  “AC”  oscillator  having  better  stability  and 
reproducibility. 


4.1.  “DC”  oscillator 

The  first  version  of  the  experimental  chaotic  oscillator,  the  “DC”  one  without  blocking 
capacitors  between  the  stages  is  presented  in  Fig.  12. 


The  Q1  -stage  is  a  buffer  and  the  Q2-stage  is  a  nonlinear  unit.  The  output  resistance  of 
the  latter  is  formed  by  the  (R1+R3)  and  R4  in  parallel.  The  linear  non-inverting  amplifier  is 
based  on  Q3  and  Q4.  We  note,  that  in  order  match  the  dc  biases  between  the  subsequent 
stages  a  p-n-p  transistor  is  exploited  for  Q3,  while  all  the  other  transistors  are  n-p-n  type  ones. 
The  delay  line  DL  (a  coaxial  transmission  cable  of  length  L)  is  coupled  to  the  matched 
resistor  R2.  The  CO  in  series  with  the  external  load  RL  is  a  blocking  capacitor,  while  CO  in 
parallel  with  the  Zener  diode  D  is  a  filtering  one.  The  diode  D,  the  auxiliary  transistor  Q5,  the 
voltage  source  V2,  and  the  resistors  RO  and  R6  serve  to  set  appropriate  dc  biases  of  the 
transistors.  The  LO  is  a  choke  used  to  block  ac  signals. 

The  specific  circuit  parameters  are  the  following:  Ro=Ri=R2=Rs=5l  Cl,  /C=/?4=75  Cl, 
R5=  11  Cl,  Re=Rg=l6  Cl,  R7= 33  Cl,  C0= 47  nF,  Ci=8  pF,  L0=4  jiH,  L=40  cm  (the  RG58U  type 
coaxial  cable;  time  delay  is  ~  5  ns/m;  Tc\ei~  2  ns),  Vi=8  V,  V2=4.5  to  5.5  V.  The  Ql,  Ql,  Q2, 
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Q4,  and  Q5  are  the  BFG520  n-p-n  type  transistors  (fti,=9  GHz),  Q3  is  the  BFR194  p-n-p  one 
(fth= 5  GHz).  The  D  is  the  BZV55-C6V2  or  similar  type  Zener  device.  In  the  RC  filter 
R=R4\\(Ri+Rl)=  47  G,  thus  the  time  constant  RC\  ~  0.4  ns.  Taking  into  account  the  time  lag  in 
the  circuit  itself  the  total  delay  time  Tdei  ~  2.3  ns,  and  the  dimensionless  delay  parameter 
T  —T de\!RC\  ~  6.  Power  spectra  of  this  “DC”  hardware  prototype  are  presented  in  Fig.  13. 

The  “DC”  oscillators  when  coupled  in  a  ring  structure  exhibit  noticeably  smoother 
power  spectrum  than  a  single  (one  module)  oscillator  (Fig.  13b  compared  with  Fig.  13a). 


Fig.  13.  Experimental  power  spectra  from  the  “DC”  chaotic  oscillator.  Vi  =8  V  and  l/2=5  V. 

(a)  single  module,  (b)  two  modules  in  a  ring  (L1+L2=L=40  cm,  C2=Ci)- 
Spectral  resolution  A/=120  kHz.  The  level  of  OdB  corresponds  to  1  pV  across  50  G  load. 


4.2.  “AC”  oscillator 

The  circuitry  of  the  “AC”  oscillator  containing  blocking  capacitors  between  the  stages 
is  shown  in  Fig.  14. 


Fig.  14.  Circuit  diagram  of  the  “AC”  chaotic  oscillator  (one  module). 


The  emitter  follower  (Q1  -stage)  is  a  buffer.  Low  dc  biased  (2  to  3  mA)  Q2-stage  plays 
the  role  of  a  nonlinear  unit.  The  output  resistance  of  the  nonlinear  unit  R5||R6  together  with 
the  capacitor  Cl  compose  the  low-pass  RC  filter.  Two-stage  circuit  (Q3-Q4)  is  a  non¬ 
inverting  linear  amplifier.  The  delay  line  DL  is  simply  a  length  L  segment  of  a  coaxial 
transmission  cable.  It  is  coupled  to  the  matched  resistors  R1  and  R15.  The  CO  are  the  blocking 
capacitors.  The  RL  is  an  external  load,  e.g.  an  input  of  a  measuring  device,  power  amplifier, 
etc.  Due  to  low  output  impedance  of  the  buffer  the  oscillator  is  not  sensitive  to  RL. 
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In  a  specific  circuit  the  parameters  are  as  follows:  R\=R5=Re=Rn=Ri5=5l  Cl,  R 2=2.1  kfi>, 
i?3=4.7  kfi,  7?4=150  Cl,  R 7=680  Cl,  l*?8=10  Cl,  Rg=R\2,=4.1  kfi,  i?io=Ri4=1.5  kf2,  Ri2=Ri6=16  Cl, 
Co=47  nF,  Ci=12  pF,  L=40  cm  (RG58U  type  coaxial  cable;  time  delay  is  ~  5  ns/m; 
Tciei~  2  ns),  Vi=5  to  6  V.  The  transistors  Q1  to  Q4  are  the  BFG520  n-p-n  type  silicon  devices 
(fth= 9  GHz).  In  the  RC  filter  R=R5\\R6=25  Cl,  thus  the  time  constant  RCi~0.3  ns  (the  parasitic 
collector-base  capacitance  of  Q1  is  small  compared  with  Ci).  Taking  into  account  time  lag  in 
the  circuit  itself  the  total  delay  time  Tdei  ~  2.3  ns,  and  the  dimensionless  delay  parameter 
T=Tdei/RC  i  ~  8.  Power  spectra  of  the  hardware  prototype  are  presented  in  Fig.  15. 


Fig.  15.  Experimental  power  spectra  from  the  “AC”  oscillator.  V^=5A  V 
(a)  single  module,  (b)  two  modules  in  a  ring  (L1+L2=L=40  cm,  C2=Ci) 

Spectral  resolution  A/=120  kHz.  The  level  of  OdB  corresponds  to  1  pV  across  50  Cl  load. 


One  can  see  in  Fig.  15  that  the  power  spectrum  in  the  case  of  two  coupled  modules  is 
much  smoother  than  in  the  case  of  a  single  module  (Fig.  15b  compared  with  Fig.  15a).  This 
experimental  observation  is  in  a  good  agreement  with  numerical  prediction  (Fig.l  la,b). 

Meanwhile  there  is  no  much  difference  in  the  power  spectra  of  the  “AC”  and  the  “DC” 
oscillators  in  the  single/two  module  cases  (see  Fig.l5a  and  Fig. 13a,  Fig.l5b  and  Fig.l3b). 


4.3.  Improved  “AC”  oscillator 


The  circuit  of  the  improved  “AC”  oscillator  stages  is  presented  in  Fig.  16. 


Fig.  16.  Circuit  diagram  of  the  improved  “AC”  chaotic  oscillator  (one  module). 
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The  role  of  the  stages  based  on  Ql,  Q2,  and  Q3-Q4  in  Fig.  16  is  the  same  as  in  the 
previous  “AC”  circuit  shown  in  Fig.  14.  The  main  circuitry  differences  are  the  following: 

(1)  The  buffer  (Ql -stage)  and  the  nonlinear  element  (Q2-stage)  are  separated  by  means  of 
the  blocking  capacitor  CO.  The  Q2  has  its  own  dc  bias  voltage  divider  R4/R5.  Though 
this  modification  requires  some  additional  circuit  elements  it  allows  one  to  set  the 
operating  point  of  the  nonlinear  element  independently  of  the  buffer  settings. 

(2)  In  the  amplifier  (Q3-Q4-stages)  large  resistors  R14  and  R18  are  inserted  in  the  emitter 
loops  to  stabilize  the  dc  conditions  of  the  transistors.  Meanwhile  the  ac  gain  of  the 
stages  depend  mainly  on  of  small  resistors  R15  and  R19  (strictly  speaking  on  the 
parallel  chains  R14||R15  and  R18||R19). 

(3)  In  contract  to  the  circuit  in  Fig.  14  the  output  stage  of  the  amplifier  (Q4-stage)  in 
Fig.  16  is  loaded  directly  with  the  DL,  the  other  end  of  which  is  loaded  with  a  matching 
resistor  R12.  This  solution  enables  one  to  increase  the  gain  of  the  stage  by  a  factor  of  2 
in  comparison  with  the  amplifier  in  Fig.  14  or  to  have  the  same  gain  but  with  a  twice 
larger  resistor  R18||R19  in  the  emitter  loop  of  the  transistor  Q4.  The  larger  values  of 
R18||R19  do  stabilize  the  gain  value  of  the  stage. 

In  a  specific  oscillator  the  circuit  parameters  are  the  following:  Ri=/?io=/?i6=6.8  k£2, 
Ri=R5=R u=Ri 7=2 .7  k£2,  Rt,=R$=R\4—R\8=100  £2,  i?4=18  k£2,  R()=R-i=R \ i~R 1 3=5  I  £2,  /?9=  1 1  £2, 
Ri5=i?i9=33  £2,  Co=47  nF,  Ci=12  pF  to  5.6  pF,  L=40  cm  to  10  cm  (RG58U  type  coaxial  cable, 
time  delay  is  ~  5  ns/m,  Tdei~  2  to  0.5  ns),  V)=8  to  9  V.  The  transistors  Ql  to  Q4  are  the 
BFG520  n-p-n  type  silicon  devices  (fh=9  GHz).  In  the  RC  filter  R=R(\\R-j=25  £2,  thus  the  time 
constant  RCi~0.3ns  to  0.14  ns  (the  parasitic  collector-base  capacitance  of  Q2  is  small 
compared  with  C\).  The  time  lag  in  the  circuit  itself  is  about  0.3  ns,  thus  the  total  delay  time 
T dd  ~  2.3  ns  to  0.8  ns  for  L=40  cm  and  L=10cm,  respectively.  Correspondingly,  the 
dimensionless  delay  parameter  x=Tdei/RC\  ~  8  at  {Ci=12  pF,  L=40  cm),  and  x=Tdei!RC\  ~  6  at 
{Ci=5.6  pF,  L=10  cm). 

Power  spectra  of  this  hardware  prototype  for  different  values  of  the  capacitances  C\,  Ci 
and  different  lengths  L  of  the  delay  line  are  presented  in  Fig.  17  and  Fig.  18. 


Fig.  17.  Experimental  power  spectra  from  the  improved  oscillator. 

(a)  single  module  (Ct=12  pF,  L=40  cm,  \/,=8.3  V), 

(b)  two  modules  in  a  ring  (L1+L2=L=40  cm,  C2=CU  l/^S.8  V) 

Spectral  resolution  A/=120  kHz.  The  level  of  0  dB  corresponds  to  1  pV  across  50  £2  load. 

One  can  see  here  again  that  the  power  spectrum  from  two-module  oscillator  (Fig.  17b)  is 
considerably  smoother  than  that  from  the  single-module  oscillator.  The  unevenness  of  the 
spectrum  is  about  10  dB  in  the  frequency  range  30  to  400  MHz  that  covers  the  VHF  band. 


S,  dB 
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Fig.  18.  Experimental  power  spectra  from  the  improved  oscillator. 

(a)  single  module  (C^=5.6  pF,  L=10  cm,  \/,=8.8  V), 

(b)  two  modules  in  a  ring  (L1+L2=L=10  cm,  C2=Ci,  \/i=8.8  V) 

Spectral  resolution  A/=120  kHz.  The  level  of  0  dB  corresponds  to  1  pV  across  50  Cl  load. 


The  two-module  oscillator  with  Ci=C2=5.6pF  and  L=  1 0  cm  exhibits  rather  smooth 
power  spectrum  (Fig.  18b).  The  unevenness  of  the  spectral  density  is  about  10  dB  in  the  UHF 
(300  to  1000  MHz)  range. 
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5.  CONCLUSIONS 

1.  The  developed  delay  line  based  oscillators  at  the  time  delays  t>  1.5 
generate  chaotic  signals.  At  larger  delays,  T  >  3  they  exhibit  hyper-chaotic 
behavior  characterized  by  multiple  positive  Lyapunov  exponents. 

2.  Two  oscillators  coupled  in  a  ring  (two-module  oscillators)  have  larger 
positive  Lyapunov  exponents,  higher  Kolmogorov  entropy,  more 
complicated  two-dimensional  structure  of  the  Poincare  sections  and 
provide  essentially  smoother  power  spectra  than  a  single  oscillator. 

3.  The  hardware  prototypes  generate  powerful  broadband  chaotic  oscillations 
in  the  very  high  frequency  (VHF:  30  to  300  MHz)  and  the  ultrahigh 
frequency  (UHF:  300  to  1000  MHz)  ranges.  The  unevenness  of  the  spectral 
density  is  about  10  dB  in  the  case  of  two-module  oscillators. 

4.  Though  the  “DC”  circuit  has  larger  bandwidths  of  the  units  there  is  no 
essential  difference  between  the  results  of  the  “DC”  and  the  “AC”  chaotic 
oscillators.  The  “DC”  version  contains  less  circuit  elements  than  the  “AC” 
one.  However,  in  the  “DC”  oscillator  it  is  much  more  difficult  to  set  the 
appropriate  dc  biases  and  to  couple  neighboring  dc  units.  In  addition,  the 
“DC”  modification  uses  the  p-n-p  type  silicon  transistors  that  commonly 
have  lower  threshold  frequencies  fth  than  the  n-p-n  ones.  Consequently,  the 
“AC”  chaotic  oscillators  have  an  advantage  over  the  “DC”  ones. 

5.  Further  research  could  be  focused  on  the  possibilities  to  increase  the 
fundamental  frequency  of  chaotic  and  hyper-chaotic  oscillations,  that  is  to 
move  or  extend  the  power  spectra  to  the  S-band  (1  to  2  GHz)  of  the 
microwave  range.  The  solutions  could  involve  different  modifications  of 
both  the  delay  line  oscillator  and  the  Colpitts  oscillator,  individual  ones  as 
well  as  the  oscillators  coupled  in  arrays. 
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